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ON  THE  SPATIAL  POWER  SPECTRUM  OF  THE  EX  B  GRADIENT  DRIFT 
INSTABILITY  IN  IONOSPHERIC  PLASMA  CLOUDS 

Experimental  studies  of  plasma  clouds  in  the  ionosphere 
[  Rosenberg,  1971;  Davis  et  al . ,  1974;  Baker  and  Ulwick,  1978] 
have  yielded  much  data  describing  ambient  ionospheric  conditions, 
e.g.,  electric  and  magnetic  fields.  The  characteristic  initial 
steepening,  elongation,  and  striation  of  E  x  B  drifting  plasma 
clouds  have  been  studied  by  applying  the  linear  theory  of  the 
E  x  B  gradient  drift  instabiltiy,  originally  developed  for  lab¬ 
oratory  gas  discharges  f  Simon ,  1963]  to  plasma  clouds  geometries 

[Haerendel  et  al.,  1967;  Linson  and  Workman,  1970;  Volk  and 
Haer endel ,  1971;  Perkins  et  al . .  1973]  .  More  recently 
Chaturvedi  and  Ossakow  [l979]  has  presented  arguments  based  on 
the  nonlinear  two-dimensional  coherent  mode  coupling  of  two 
Fourier  modes  of  plasma  cloud  density,  following  the  work  of 
Rognlien  and  Weinstock  [  1975  ],  to  explain  the  nonlinear 
stabilization  and  resultant  saturated  amplitudes  and  power 
spectra  of  the  long  wavelength  E  x  B  gradient  drift  instability 
in  ionospheric  plasma  clouds. 

Numerical  simulations  \ Zabusky  et  al . ,  1973;  Lloyd  and 

Haerendel ,  1973;  Goldman  et  al . ,  1974;  Doles  et  al . ,  1976; 
Ossakow  et  al . ,  1975,  1977]  of  the  interaction  of  barium  plasma 

clouds  with  the  ionosphere  have  reproduced  not  only  many  of  the 
gross  observational  features  of  plasma  cloud  evolution,  but 
also  their  spatial  power  spectra  fScannapieco  et  al.,  1976], 
minimum  striation  scale  sizes  [McDonald  et  al . .  1978,  1980] 
and  outer  scale  size  or  correlation  length  [ Keskinen  et  al  .  , 

1980a ;  Keskinen  and  Ossakow,  1 980  ] .  In  addition  numerical 
Manuacript  aubmitted  February  3, 1981. 
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simulation  studies  of  the  E  x  B  gradient  drift  instability  in 
local  unstable  regions  of  ionospheric  plasma  clouds  [ Keskinen 
et  al. ,  1980b]  have  yielded  spatial  power  spectra  and  satur¬ 
ated  wave  amplitudes  consistent  with  experimental  values. 

However,  to  date,  an  analytical  description  of  the  two- 
dimensional  spatial  power  spectra  of  the  E  x  B  gradient  drift 
instability  in  ionospheric  plasma  clouds  over  a  wide  range  of 
wavelengths  has  not  been  discussed  in  detail.  In  this  paper 
we  show  that  the  spatial  power  spectra  of  ionospheric  plasma 
clouds  recently  computed  from  numerical  simulation  studies  and 
experimental  observations  are  consistent  with  conservation  laws 
implied  by  the  fundamental  fluid  equations  modeling  the  E  x  B 
gradient  drift  instability  in  ionospheric  plasma  clouds. 

For  wavelengths  much  greater  than  the  ion  gyroradius 
(approximately  10  meters  for  Ba+  ions  in  the  twilight  F  region), 
the  dynamics  of  the  plasma  cloud  and  background  ionosphere  can 
be  studied  in  the  fluid  approximation  fvb'lk  and  Haerendel, 

1971;  Perkins  et  al .  ,  19  73  ;  Zabusky  et  al .  ,  1973;  Ossakow  et  al., 
1975].  For  large  clouds  (large  magnetic  field  line  integrated 
Pedersen  conductivity  compared  with  that  of  the  background 
ionosphere),  the  cloud  interaction  with  the  background  iono¬ 
sphere  (second  level)  can  be  neglected  f Haerendel  et  al.,  1967]. 

Furthermore,  due  to  the  very  high  conductivity  along  the  mag- 

_  6 
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field  line  integrated  (two-dimensional)  model  is  justified 
r Volk  and  Haerendel.  1971;  Perkins  et  al . ,  1973].  In  reality, 
an  artificially  injected  plasma  cloud  will,  initially,  be  two- 
dimensional  in  the  plane  perpendicular  to  the  magnetic  field. 

But  a  two-dimensional  plasma  cloud  of  finite  extent  is  not 
stationary  f Dungey ,  1958]  and  will  distort  as  it  convects.  How¬ 
ever,  for  striation  scale  sizes  less  than  the  ambient  cloud  Pe¬ 
dersen  conductivity  gradient  scale  length  and  on  the  faster  time 
scale  on  which  striations  develop,  the  assumption  of  one-dimen¬ 
sionality  is  justified.  As  a  result,  in  what  follows  we  confine 
our  attention  to  local  regions  of  initially  slablike  (one¬ 
dimensional)  ionospheric  plasma  clouds.  This  will  also  facili¬ 
tate  comparlsion  with  numerical  simulation  studies. 

By  adopting  a  Cartesian  coordinate  system  (x,y,z)  with 
magnetic  field  Bz,  ambient  electric  field  EQy,  and  after  trans¬ 
forming  to  a  frame  drifting  with  velocity  Vq  =  (cEq/B)x,  the 
two-dimensional  model  equations  for  the  magnetic  field  line  in¬ 
tegrated  plasma  cloud  Pedersen  conductivity  £(x,y)  and  the  self- 
consistent  plasma  cloud  electrostatic  potential  Sq-(x,y)  can  be 
written 


+  -  z  x  V<$c h-VE  -  DV  ^  £ 

*J  t  D 


V  •  (  E  V  5tp)  =  Eq  •  VZ 


where  V <5 v  =  ~  E(x,y)  +  Eq  ,  E(x,y)  the  total  electrostatic  field, 
c  is  the  speed  of  light,  V  =  (3/3x,  9/3y)  and  D  is  the  cross¬ 
field  diffusion  coefficient  [Perkins  et  al.,  1973]  given  by 
2 ( vg /ftg ) (c kgT/ eB )  with  T  the  ion  and  electron  temperature. 
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vg  the  sum  of  electron  collision  frequencies  with  plasma  cloud 

ions  and  ambient  neutrals,  k„  is  Boltzmann's  constant,  and  Q 

B  ’  e 

the  electron  gyro f requency .  All  other  symbols  retain  their 
conventional  meaning.  Equation  (1)  results  from  the  magnetic 
field  line  integration  (along  the  z-direction)  of  the  ion  con¬ 
tinuity  equation,  while  equation  (2)  is  derived  from  current 
conservation  V*J_  =  0. 

By  linearizing  (1)  and  (2)  and  assuming  fluctuations  in 
magnetic  field  line  integrated-Peder sen  conductivity  and  cloud 
potential  6  Z ,  6cp  exp  ^i(kyy  +  k^x)  +  Y^tJ  w^t^1  k*B  =  0,  kL>>l, 
one  finds  the  usual  J2  x  15  growth  rate 

Yk  =  <cEo/BL)(ky/k)2  -  Dk2  (3) 

where  k2  =  k  2  +  k  2 ,  L  =  31nE  /3x.  For  cE  / B  =•  100  m/sec, 
x  y  *  o  o  *  » 

L  =  6  km,  D  =  lm  /sec,  the  critical  wavelength  \c  (y^  *0)  =60  m. 

We  will  now  show  that  several  features  of  the  spatial  power 
spectra  of  the  E  x  B  gradient  drift  instabiltiy  in  ionospheric 
plasma  clouds  as  observed  in  experiments  and  recent  numerical 
simulations  can  be  derived  by  considering  the  conservation  laws 
implied  by  (1)  and  (2).  We  can  write  (1)  and  (2),  after  making 
the  separation  E(x,y)  =  Eq(x)  +  6E(x,y),  as  follows: 

+  §  z  x  V  69 • V  E  +  §  z  x  V  69 • V  6  E  -  DV  2  6  E  (4) 
at  d  —  OB 

EqV2<59  +  6  E  V  2  69  +  V  E  o  •  V  69  +  V  6  E  •  V  69  =  Eq«V6E  (5) 

We  first  multiply  (4)  by  <5  E  and  integrate  over  all  x  and  y 
finding 
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Ihf* 


d2x  (5E)2  + 


if  <2* 


6 E  z  x*?  6cp  •  VE 


+  ^ J“ d2x  6E  z^  xV<$(p 


•75E  -  nj< 


d2x  <$EV2<5E 


(6) 


where  d2x  =dxdy.  The  third  term  on  the  left  hand  side  vanishes 
s  ince 


J" d2  x  5E  £xV6o*V6E 

if2'  7 

=  Id  x  _z  x  V6(p*V(6E)2 

=  f/d2x  V^(dE)2zxV64<] 


a  0  if  (6£>— ►  ()  as  x , y- 


since  V*|  (6E)2  z  xV6cp 


]**  X 


V  •  V  (  6  E  )  2  .  After  making  the  follow¬ 


ing  Fourier  expansion  of  £E  and  6q< 


6E (x, y) 
6^(x,y) 


d2k’  6E,  ' 
k 

d2k  6cpk 


ik’  *x 

e  —  - 

ik  ■  x 

e  —  — 


the  second  term  on  the  left  hand  side  of  (6)  can  be  written 


|  f  d2x J" d2kr  6ER’  e  i-'  zx7  (  /d 2  k  6(pu  e  ]*VE 


0 


-  j^fd2kdzk' 


i  z  x  k  •  V  £  6  E 

—  —  o 


-) 

k'  5*  A2x  «  ICJt- 


d2k  d2k’  i  X  x  k*  72  52.’  dv.5(k’+k) 


o  k 


£  f  d2k  i  z  x  k-VE  5 E  ,  6q,. 

B  J  ~  o  -k  '♦'k 


where  6 (x)  is  the  two-dimensional  delta  function.  The  first 
and  fourth  terms  in  equation  (6)  are  reduced  in  a  similar 
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manner 


and 


1  Tif  '*2»<8E>2  -  5  Tif  "2k  l8£kl 


d2  x  5 1  V26E 


i 


d2k  k2l«kl2 


All  that  remains  is  to  find  an  expression  for  6<f>k  to 
evaluate  the  second  term  on  the  left  hand  side.  Since  we  are 
considering  only  quadratic  nonlinearities  only  a  linear  relation 
of  <5<pk  to  is  needed  and  one  can  use  the  linearized  version 

of  equation  (S)  to  find  in  terms  of  ££k. 

I  V26<p  +  VE  *V6<p  -  £  *V6 £  (5a) 

O  o  o 

Inserting  Fourier  expansions  of  6£  and  6<p  into  (5a)  we  find 
in  a  straightforward  manner 

6vk  as  -  (ik-Eo/k2ro)62:k  (7) 


Substituting  into  the  Fourier  transformed  version  of  the  second 
term  on  the  left  hand  side  of  equation  (6)  gives 


(z  x  k*VE  k 
—  o  — 


•E  /k2E  ) 
— o  o 


Equation  (6)  can  than  be  written  as 


_1 

2 


'/ dak  Tk  pk 


(S) 


where  P.  =  J  6  E  | 2  -  |6E(k  ,k  )|2  and  y.  is  the  linear  growth 
k  x  y  x 

rate  given  by 
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Yv  =  %  (k  •  E  z  xk • V£  /k2Z  )  -  Dk2 
JE  B  —  — o  —  —  o  o 

=  (cE  /BL)  cos  20  -  Dk2 
o 

where  L  1  -  ( 1 / Z  ) ( D Z  /3x)  and  0  is  the  angle  defined  by  k  and 
o  o  — 

E  .  We  note  that  the  nonlinear  term  does  not  enter  into  (8); 

— o 

the  nonlinear  term  thus  conserves  energy  and  may  be  thought  of 
as  a  bridge  between  the  stable  and  unstable  regions  of  wave- 
number  space.  This  property  is  also  shared  by  several  other 
nonlinear  plasma  physics  problems,  e.g.,  the  two-dimensional 
electrostatic  guiding  center  plasma  f Montgomery ,  19763,  ion- 
acoustic  shock  formation  [Ott  et  al . ,  1973],  and  two-dimensional 
turbulence  in  the  equatorial  electrojet  f Keskinen  et  al.,  1979; 
Sudan  and  Keskinen,  1977  ,  1 9  7  9  3  - 

In  the  steady  state  (3/3t  =  0)  we  can  write  (8)  in  polar 
coordinates  (kz  ®kx2  +  ky2  »  tan  ©  =  kx/ky)  assuming  the 
spatial  power  spectrum  extends  from  kmin  to  ^max 


.  k 

J  k  J  J  o 

min 


2  IT 


d0  Pk 


/k  2u 

max  r 

dk  k 

Jo 


d0  *k  Pk 


(9) 


with  y®  =  (cE  /BL)  cos20  and  y,  =  Dk2 .  As  a  result  the  total 
k  o  k 

power  generated  in  the  unstable  range  of  w-.venumbers  between 

k  .  and  k  ,  the  critical  wavenumber,  must  be  balanced  by  the 
min  c 

power  dissipated  in  the  stable  regime  between  k£  and  ^max  under 
steady  state  conditions. 


Our  approacn  will  be  to  make  an  ansatz  for  P(k,0)  and  show 
its  consistency  with  (9).  Previous  numerical  simulation  studies 


f  Scannapieco  et  al. .  1976;  Keskinen  et  al. .  1980a;  Keskinen  and 
Ossakow ,  1980j  have  shown  that  the  one-dimensional  transverse 
integrated  spatial  power  spectra  of  magnetic  field  line  inte¬ 
grated  Pedersen  conductivity  perturbations  in  the  plane  per¬ 
pendicular  to  the  magnetic  field  of  initially  slablike  iono¬ 
spheric  plasma  clouds  can  be  represented  by  a  power  law,  i.e., 

/dk  |  6  E  (k  ,  k  )|2  =  (k  2  +  k  2)_ny/2 

x  1  x  y  '  y  oy 

and 

/  <*ky  l^<>VVl2  '  <kx2  +  ko^)"”k'  2 

with  n  =  2,  n  -  2-3,  and  2ir/k  the  outer  scale  size  where  the 
x  y  o 

x-axls  denotes  the  direction  of  plasma  cloud  E  x  B  drift.  These 
power  spectral  Indices  are  in  agreement  with  those  derived  from 
experimental  observations  [Baker  and  Ulwick.  1978;  Kelley  et  al., 
1979]  of  ionospheric  plasma  clouds.  Recently,  the  angular  de¬ 
pendence  of  P(k,0)  has  been  studied  f Keskinen  and  Ossakow,  198C] 
in  the  nonlinear  steady  state  regime  for  initially  slablike 
ionospheric  plasma  clouds  under  different  sets  of  initial  con¬ 
ditions.  In  polar  coordinates,  these  studies  show  that  P(k,0) 
is  anisotropic  with  maximum  power  along  and  near  the  y-direction 
(perpendicular  to  Eq  x  B  direction)  and  suggest  that 

P (k,  Q  )  £:  P  cos  m0  (1  +  (k/k  )2)-(n+l)/2 
o  o 

with  m>0,  P  s  const,  and  2rr/k  is  the  outer  scale  size.  We 
o  o 

note  that  this  form  for  P(k,0)  maximizes  in  the  linearly  most 
unstable  direction  (0  -  0  where  0  is  defined  by  k  and  Eq).  This 
anisotropy  is  a  reflection  of  the  fact  that  nonlinear  plasma 
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cloud  striations  are  elongated  much  more  in  the  Eq  x  B 
direction  and  this  result  is  also  observed  in  local  simulations 
of  ionospheric  plasma  clouds  f Keskinen  et  al.,  1980b].  The 
anisotropic  nature  of  the  two-dimensional  spatial  power  spectra 
of  ionospheric  plasma  clouds  is  also  consistent  with  arguments 
put  forth  by  Chaturvedi  and  Qssakow  [l 97  9]  in  discussing  the 
nonlinear  stabilization  of  the  E  x  B  gradient-drift  instability 
in  ionospheric  plasma  clouds.  In  addition,  this  anisotropy  is 
also  observed  in  experimental  plasma  cloud  studies  [M.C.  Kelley, 
private  communication,  198c]. 


We  substitute  this  ansatz  for  P(k,0)  into  the  left  hand 


side  of  (9)  obtaining 


WBL)/ 

J  O 


d0  cosm+2  of  dk  k  ( 1  +  ( k/k  )2)_(n+D/2 


and  we  have  assumed  k>k  ,  k  <  k  <  k  and  n>  1 .  Similiarly  the 

cominoc 

right  hand  side  of  (9)  gives 
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4 


P  D  2/7i 
o 


r<f+l) 


k  3-n 

(_jn_aX) 


We  note  that  since  the  left  and  right  hand  sides  of  (9) 
must  be  positive,  the  spectral  index  l<n<3.  This  implies  that 
the  steady  state  spatial  power  spectrum  can  be  represented  by 
a  power  law.  Equating  the  power  generated  in  the  unstable  re¬ 
gion  to  the  power  dissipated  and  solving  for  the  outer  scale 

wavenumber  k  we  find 
o 


3  -n 
n-  1 


r<!+2) 

r  <;+l )  r  (  S±2) 


,cE  ,BL 

(  o/  ) 

Dk 

max 


(10) 


For  an  inverse  power  law  with  spectral  index  n  =  2  and  a  cos^G 


angular  (m  =  2)  dependence  we  find  kQ*L 


This  scaling  of  out¬ 


er  scale  wavenumber  k  with  the  initial  Pedersen  conductivity 

o 

gradient  scale  length  when  na2-3  is  consistent  with  recent  numeri¬ 
cal  simulation  studies  of  the  E  x  B  gradient  drift  instability 
in  ionospheric  plasma  clouds  [ Keskinen  et  al .  ,  1980a;  Ke sk in en 

and  Ossakow,  1980].  Furthermore,  by  taking  D  =  1  m2/sec, 

L  =  6  km,  V  =  cE  /B  =  100  m/sec,  2i\/k  ~  20  m  f McDonald 

et  al . ,  1978;  Baker  and  Uiwick,  1978]  we  find  from  (10)  that 

2ir/k  =200  m. 
o 

We  further  note  that  these  results  may  also  be  applicable 
to  the  natural  ionosphere  since  if  one  replaces  the  terms 
DVz8Z  and  Eq*V6Z  in  (4)  and  (5)  with  -v^6I  and  v  V  6  E  f 
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respectively,  where  v  is  a  recombination  rate,  £  denotes 
gravity,  and  the  collision  frequency  with  neutrals,  then 

(4)  and  (5)  are  similar  to  the  equations  that  describe  the  col- 
lisional  Ray leigh-Taylor  instability  [Ossakow  et  al .  ,  1979]  in 
the  equatorial  F  layer.  The  adaptation  of  this  method  to  the 
coliisional  Rayleigh-Taylor  instability  and  current  convective 
instability  is  now  in  progress. 
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OICY  AT’N  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSI’Y  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  308 
LIVERMORE,  CA  99550 

OICY  A’TN  doc  CON  for  tech  INFO  oept 
OICY  A”N  DOC  CON  FOR  --589  ft.  0” 
OICY  A"N  OOC  CON  FOR  l-31  ft.  HAGEft 
OICY  ATTN  JOC  CON  FOR  l-96  F.  SEWARD 


LOS  alamos  SCIENTIFIC  laboratory 
P.  0.  BOX  1663 


alamos, 

, 

87595 

oicy 

at-* 

DOC 

CON 

FOR 

J.  WOLCOTT 

OICY 

attn 

DOC 

CON 

FOR 

R.  F.  TASCHEK 

OICY 

a 

DOC 

CON 

FOP 

E.  JONES 

OICY 

AT"TN 

DOC 

CON 

FOR 

J.  MALIK 

OICY 

ATTN 

DOC 

CON 

FOR 

R.  JEFFRIES 

OICY 

A~N 

DOC 

CON 

FOR 

J.  ZI** 

OICY 

at-* 

DOC 

CON 

FOR 

P.  KEATON 

OICY 

ooc 

CON 

FOR 

D.  wES'ERi®.’ 

SATE)  1 A  LABORATORIES 
P.  0.  BOX  5100 
ALBUQUERQUE,  KM  87115 

OICY  ATTN  DOC  CON  FOR  J.  MARTIN 
OICY  at-*n  DOC  CON  for  w.  BROWN 
oicr  a”n  doc  con  for  a.  tvornbrough 
OICY  A”n  doc  CON  for  T.  WRIGM’ 

OICY  attn  doc  CON  FOR  D.  OAFLGREN 

OICY  a”n  doc  con  for  3191 

OICY  AT-H*  DOC  CON  fqr  SPACE  PROuEC’ 


sandia  -abora’ories 

LIVERMORE  laboratory 
p.  0.  BOX  969 
LIVERMORE,  CA  99550 

OICY  A’-N  OOC  con  for  B.  MiRPfEY 
OICY  Ar*N  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MIL  I 'ARY  APPLICATION 
DEPARTIN’  of  f>«R&y 
WASMING’ON,  D.C.  20595 

OICY  a”n  DOC  CON  'OR  0.  GALE 


OTHER  iOvERKMEN' 


SAMSO/SK 

P.  0.  BOX  92960 

WORLDWAY  POS’AL  CEN’ER 

LOS  ANGELES,  CA  90009 

OICY  A”N  S«A  (SPACE  COMM  SYS’EMS)  M.  CLAVIN 


central  INTELLIGENCE  AGENCY 
attn  RO/SI,  RM  5G9g.  HQ  BLOG 
MASMING’ON,  D.C.  20505 

OICY  A”N  OSI/PS1D  RM  5F  19 


5APISO/FM 

NOR’ON  APB,  CA  92909 

(MINUTEMAN) 

OICY  AY’S  tPA.  L’C  KEKMtDY 


DEPARTIN’  OF  COMMERCE 
NA’IONAL  BUREAU  OF  j'AKOAROS 
wASHlNG’ON,  D.C.  20239 

(ALL  CORRES.  A"n  SEC  OFFICER  FOR) 
OICY  a”n  «.  MOORE 


DlV 
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ins'iyu'e  for  telecom  sciences 
national  reucowutiCAriONS  i  i**o  *omin 

HU.OU,  CO  10 SOS 

OICY  A'TN  A.  JEAN  (UNCLASS  0M.Y) 

OICY  attn  w. 

oicy  attn  o.  crohbie 

OICY  ATTN  t  .  KMT 

NATIONAL  OCEANIC  i  ATMOSPHERIC  AOMIN 
ENVIRONMENT*.  RESEARCH  LABORATORIES 

oepartknt  op  compcrci 
boulder,  co  80)02 

OICT  attn  a.  MUM 

Oicr  ATTN  AERONOMY  L A«  G.  REID 


OEPARTKWT  OP  OEFENSE  C  ONT  RAC ’ORS 


AEROSPACE  CORPORA' ION 

P.  0.  BOX  929S7 

LOS  AMGELES,  CA  90009 


OICT 

attn  i. 

•arpumcel 

OICT 

Arm  r. 

SALMI 

OICT 

AT’N  V. 

vOSEPhSON 

OICT 

A'tn  s. 

BONER 

OICT 

AT'N  N. 

stocknell 

OICT 

AT’N  0. 

OLSEN 

OICT 

attn  SPPA  FOR  pnw 

analytical  SYSTEMS  enciphering  COM 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  0180) 

OICY  attn  RADIO  SCIENCES 


BERKELEY  RESEARCH  ASSOCIATES,  INC. 

P.  0.  BOX  91) 

BERKELEY,  CA  94  701 

OICY  attn  j.  WORKMAN 

BOEING  COMANY,  TK 
p.  0.  BOX  5707 
Seattle.  wA  98124 

OICY  A'tn  g.  CEIS'ER 
OICY  ATtn  0.  -AJRRAY 
OICY  A’»N  G.  MALL 
OICY  attn  j.  XEPPHY 

CALIFORNIA  AT  JAN  OIEGO,  JNIV  OP 

P.O.  Box  6049 
San  Diego,  CA  92106 

BRCMN  ENGINEERING  COM  ANY,  INC. 
COMINGS  RESEARCH  PARK 
HUATSvlLLE,  al  35807 

OICY  attn  ROMO  A.  oetlMRIS 

Charles  s’ark  draper  laboratory,  Inc. 
555  TEO*OL06Y  SQUARE 
CAMRIOGE,  MA  021)9 

OICY  ATTN  o.  8.  COX 
OICT  ATTN  j.  P.  GILMORE 


COMSAT  LABORATORIES 

l i  nth ( cum  *o*o 
CLAMSauRG,  TO  207)4 
OICT  ATTN  G.  HYDE 

CORNELL  UNIVtRSITY 

oepartmnt  op  electrical  ewimering 

I THACA,  NY  141)0 

OICY  ATTN  3.  T.  PARLEY  UR 


ElECtROSPACE  SYSTEMS,  INC. 

BOX  1)59 

RICHAROSON,  Tx  75080 

OICY  A»tn  h.  lOGStON 

OICY  attn  SECURITY  (PAUL  PHILLIPS) 

ESl  INC. 

49S  JAVA  DRIVE 
SUNNYVALE,  CA  94086 

OICY  ATTN  J.  ROBERTS 
OICT  AYYN  janes  MARSHALL 
OICY  A'TN  C.  W.  PRETTIE 

FORD  AEROSPACE  6  COMMUNICATIONS  CORP 
59)9  PARIAN  WAY 
PALO  AL'O,  CA  94)0) 

OICY  AT'N  j.  T.  MAYT|NGLEY 

general  electric  company 
space  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GOOOARO  BlvO  XING  OF  PRUSSIA 
P.  0.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICT  attn  m.  h.  BORTNER  SPACE  SCI  lab 

gephral  electric  company 

P.  0.  BOX  1122 
SYRACUSE,  NY  15201 

OICY  AT’N  F.  REISER* 

GEPHRAL  ELECTRIC  COMPANY 

TEMPO -CENTER  for  ADVANCED  S'UOIES 

816  S'A'E  street  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  95102 


OICT 

ATTN 

DAS  1  AC 

OICY 

at*n 

don  chandler 

OICY 

AT^Tn 

tom  BARRETT 

OICY 

AT-TN 

TIM  S'EPHANS 

OICY 

A^N 

WARREN  S.  KNAPP 

OICY 

A^N 

WILLIAM  MCNAMARA 

OICY 

AT^ 

B.  GAMHI.L 

OICY 

AT-*  SrANTON 

GENERAL  ELECTRIC  tech  SERVICES  CO.,  INC. 

COURT  S'REET 
STRACuSE,  NY  15201 

OICY  AT'N  G.  MILLMAN 

general  research  corporation 

SANTA  BARBARA  DIVISION 

P.  0.  BOX  6770 

SANTA  BARBARA,  CA  95111 

OICY  A'tn  JOHi  1SE  JR 
OICY  AT'N  JOEL  GARSAR1NO 


GEOPHYSICAL  INSTITUTE 
lAlveRSI'Y  OP  ALASKA 
FAIRBAMLS,  AK  99701 

(All  CLASS  Arm;  SECURITY  OFFICER) 
OICY  AT*n  T.  N.  DAVIS  OJNCL  OR.Y) 

OICY  attn  kal  SROP*  (UNCL  0M.T) 

OICT  attn  TECHNICAL  library 

GTE  STL  VANIA,  INC. 

ELECTRONICS  STS’IMS  grp-eastern  DIV 
77  A  street 
HEOHAH,  MA  02194 

OICT  attn  marshal  cross 

Illinois,  university  op 
departin’  op  electrical  enciphering 

JRSANA,  IL  6180) 

OICT  attn  x.  TEH 

Illinois,  university  op 
107  COBLE  HALL 
101  S.  WRIGHT  STREET 
jftSAM,  It  60680 

(ALL  COR RES  AT'n  SECURITY  SUPERVISOR  *OR) 
OICY  attn  x.  TEH 
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institute  fo?  osfsnsc  analyses 

400  ARMY  -NAVY  DRIVE 
ARLINGTON,  VA  22202 

01CY  ATTN  J.  M.  A£ IN 
01CY  ATTN  ERNES r  BAUER 
OICY  ATTN  HANS  WQCFHARD 
OICY  ATTN  jOEl  BENGSTON 

HSS,  INC. 

2  ALFR60  CIRCLE 
BEDFORD,  HA  01730 

01CY  Ar'rN  DONALD  HANSEN 


1NTL  TEL  6  TELEGRAPH  corporation 
500  *ASHlNGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  AT’N  technical  library 


JAYCGR 

1401  CAN  I  NO  DEL  HAR 
DEL  HAR,  CA  32014 

OICY  ATT*  s.  R.  GOLDMAN 


JOHNS  -OP<lNS  JNI'.ERSITY 
APPLIED  PHYSICS  '-A80RAT0RY 
JOHNS  OPKINS  ROAD 
LAUREL,  MD  20810 

OICY  AY’n  OOCUMENt  LIBRARIAN 
Oicr  attn  Thomas  po'ehia 
OICY  ArT*  jOhN  OASSOULAS 

LOCKHEED  MISSILES  i  SPACE  CO  INC 
P.  0.  BOX  504 
SUNNYVALE,  CA  34088 

OICY  ATtn  OEpT  60-l2 
OICY  ATtn  D.  R.  CHURCHILL 

LOCKHEED  MISSILES  AND  SPACE  CO  INC 
5251  HANOVER  S'rR££T 
PALO  AL-O,  CA  94304 

OICY  MARGIN  MALt  XPt  52-20 

OICY  attn  richARD  G.  JOHNSON  DEPT  52-12 

oicy  a^n  if^of  ogpr  52-12 

KAHAN  SCIENCES  CORP 
P.  0.  BOX  7463 
COLORADO  SPRINGS,  CO  80935 
OICY  a^n  r.  MEAGHER 

UNCASI”  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  A^tn  jRW|N  JACOBS 

M.l.T.  LINCOLN  LABORATORY 
p.  0.  BOX  73 
LEXINGTON,  ha  02173 

OICY  A*’*  qavIO  H.  TOWLE 

oicr  at^n  p.  maldron 

OICY  ATTN  l.  LOUGMLIN 
Oicr  ATTN  0.  CLARK 

HART  IN  HAR l ETtA  CORP 
ORLANOO  DIVISION 
P.  0.  BOX  5837 
ORLANOO,  CL  32805 

OICY  ATTN  R.  HEFFNER 

HCOONNELL  DOUGLAS  CORPORATION 
5301  90LSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATtn  n.  HARRIS 
OICY  ATTn  j.  moule 
OICY  attn  GEORGE  HROZ 
OICY  attn  OLSON 
OICY  Attn  R.  w.  HALPRIN 

OICY  AT-N  technical  library  services 


HISS  ION  RESEARCH  CORPORATION 
735  $*ATE  STREET 
SANTA  BARBARA,  CA  93101 
OICY  AT^  p.  FISCHER 
OICY  ArTN  tf.  F.  CREVIER 
OICY  attn  steven  L.  GUTSCHE 
OICY  attn  SAPPE^IELO 
OICY  ATtn  r.  8QGUSCM 
OICY  Attn  r.  HENDRICK 
OICY  AT’N  RALPH  <IlB 
OICY  attn  DAVE  SOWLE 
oicy  AT*N  f#  fajEN 
OICY  attn  m.  SCHE10E 
OICY  attn  COMRAD  l.  lONGMIRE 
OICY  ATTN  wARREN  a.  schlje^er 

MITRE  CORPORATION,  THE 
P.  0.  BOX  208 
BEDFORD,  HA  01730 

OICY  ATTN  JOHN  HORGANSTERN 
OICY  attn  g.  HARDING 
OICY  AT-n  C.  e.  CALLAHAN 

MItre  CORP 

WCS-GATE  RESEARCH  PARK 
1820  DOLLY  MADISON  3LVD 
MCLEAN,  VA  22101 

OICY  A^N  W.  HALL 
OICY  ATtn  w.  pOStER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELO  BLVD. 

SANTA  MONICA,  CA  90404 

OICY  At-n  E.  C.  field  JR 

PENNSYLVANIA  S^ATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 

318  electrical  engineering  east 

UNIVERSITY  PARK,  °A  1S802 

(NO  CLASSIFIED  to  ^1$  ADORESS) 

OICY  attn  ionospheric  research  LAB 

PHOTOMEtricS,  INC. 

442  HAR RETT  R0A0 
LEXINGTON,  HA  021 75 

OICY  ATtn  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  3027 
BELLEVUE,  wA  38009 

OICY  Attn  e.  j.  FREHOUW 

physical  dynamics  inc. 

P.  0.  BOX  1036/ 

OAKLAND,  CA.  94610 

ATTn:  A.  THOMSON 


R  4  D  ASSOCUTES 
0.  BOX  9695 

MARINA  DEL  RET,  CA  90291 

OICY  ATtn  FORREST  GILMORE 
OICY  ATTN  3ft  Y  AN  GABBARD 
OICY  AT^N  WILLIAM  B.  WRIGHT  jR 
OICY  ROBERT  f.  LELEVIER 

OICY  attn  william  j.  karzas 
OICY  AT^n  m.  dry 
OICY  attn  c.  MACDONALD 
OICY  ATtn  ft.  TURCO 

rano  corporation.  The 
1700  MAIN  S^REgT 
SANTA  MONICA,  CA  90406 

OICY  4TTN  CULLEN  CRAIN 
OICY  ATTN  ED  B6DR02 I AN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  EM)  AVENUE 
YORK,  NY  1002  3 
OICY  attn  VINCE  TRAPANI 


20 


SCIENCE  APPLICATIONS,  INC. 

P.  0.  BOX  2331 
LA  JOLLA,  CA  92038 

01CT  ATTN  LEWIS  M.  L1NS0N 
01CY  ATTN  DANIEL  A.  H*UN 
OICT  ATTN  0.  SAOfS 
01CT  ATTN  E.  A.  STAKE* 
OICT  ATTN  CURTIS  A.  SMITH 
OICT  ATTN  JACX  NC00U1ALL 


Raytheon  Co. 

528  Boston  Post  Road 
Sudbury,  MA  01776 

01CY  ATTN  Barbara  Adams 

Science  Applications,  Incorporated 
1710  Goodridge  Dr. 

McLean,  VA  22102 

Attn:  J.  Cockayne 


SRI  INTERNATIONAL 
533  RAVENSWOOO  AVENUE 


MEM.0  PARK, 

7  CA 

94025 

OICT 

ATTN 

DONALD  NEEL  SON 

OICT 

ATTN 

ALAN  BURNS 

OICT 

ATTN 

5.  SMITH 

OICT 

ATTN 

L.  L.  COBB 

OICT 

ATTN 

DAVID  A.  JOHNSON 

OICT 

ATTN 

WALTER  G.  CmESNUT 

OICT 

AT^N 

CHARLES  L.  RINO 

OICT 

AT-N 

WALTER  JAYE 

OICT 

ATTN 

M.  BARON 

OICT 

ATTN 

RAY  L.  LEADABRAND 

OICT 

ATTN 

G.  CARPENTER 

OICT 

ATTN 

G.  PRICE 

OICT 

ATTN 

J.  PETERSON 

OICT 

ATTN 

R.  HAKE,  JR. 

OICT 

ATTN 

V.  GONZALES 

OICT 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
73  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICT  ATTN  W.  9.  BOQUIST 


Lockheed  Missile  &  Space  Co.,  Inc. 
Huntsville  Research  &  Engr.  Ctr. 
4800  Bradford  Drive 
Huntsville,  Alabama  35807 
Attn:  Dale  V.  Davis 


TRW  DEFENSE  t  SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDOAOO  BEACH,  CA  90278 

OICT  ATTN  R.  <•  PLEBUCH 
OICT  AHN  s.  ALTSCMJLER 
OICT  ATTN  0.  DEE 

VISIDY*,  INC. 

19  THIRD  AVENUE 

NORTH  WEST  INDUSTRIAL  PARK 

BURLINGTON,  MA  01803 

OICT  ATTN  CHARLES  HUMPHREY 
OICT  ATTN  J.  W.  CARPENTER 
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STRATEGIC  TECHNOLOGY  OFFICE 
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DEPARTfCNT  OF  THE  NAVY 
WASHINGTON,  D.C.  20560 

OR.  T.  CZUBA 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  2057$ 


DR.  P.  MANGE 
OR.  R.  FCIER 

OR.  E.  SZUSZCZEW1CZ  -  C30E  6127 

OR.  TIMOTHY  COFFEY  -  CODE  4700  (2$  COPIES) 

DR.  S.  OSSAKOW  -  COOE  4780  U88  COPIES) 

OR.  J.  GOODMAN  -  COOE  7560 


SCIENCE  APPLICATIONS,  INC. 
1250  PROSPECT  PLAZA 
LA  JOLLA,  CALIFORNIA  92057 


OR.  D.  A.  HAH. IN 
OR.  L.  LINSON 
OR.  D.  SACHS 


DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL  LABORATORY 
NOAA 

BOULDER,  COLORADO  80502 

OR.  A.  GLENN  JEAN 
DR.  G.  W.  ADAMS 
OR.  D.  N.  ANDERSON 
DR.  K.  DAVIES 
OR.  R.  F.  DONNELLY 


HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MASS.  02158 


OR.  M.  B.  MCELROY 
DR.  R.  LINOZEN 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PENNSYLVANIA  16802 


DR.  J.  S.  NISBET 
OR.  P.  R.  ROHRBAUGH 
DR.  D.  E.  BARAN 
DR.  L.  A.  CARPENTER 
DR.  M.  LEE 
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UNIVERSITY  OF  CALIFORNIA,  LOS  ANGELES 

405  HILLGARD  AVENUE 

LOS  ANGELES,  CALIFORNIA  90024 

OR.  F.  V.  CORONITI 
OR.  C.  KENNEL 


A.  F.  GEOPHYSICS  LABORATORY 
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OR.  T.  J.  KENESHEA 
OR.  J.  AARONS 


OFFICE  OF  NAVAL  RESEARCH 
BOO  NORTH  QUINCY  STREET 
ARLINGTON,  VIRGINIA  22217 


UTAH  STATE  UNIVERSITY 
4TH  N.  AND  8th  STREETS 
LOGAN,  UTAH  84322 

DR.  P.  M.  BANCS 
OR.  R.  HARRIS 
OR.  V.  PETERSON 
OR.  R.  PEG ILL 
DR.  K.  BAKER 
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CORNELL  UNIVERSITY 
ITHACA,  NEW  YORK  14850 


coftwccrc 

naval  ELECTRONICS  LABORATORY  CENTER 
SAN  01  EGO,  CALIFORNIA  92152 

DR.  M.  BLEIWEISS 
OR.  I.  ROTHM.LER 
OR.  V.  HILOEBRAM) 

MR.  R.  ROSE 


OR.  W.  E.  SWARTZ 
OR.  R.  SUDAN 
DR.  0.  FARLEY 
DR.  M.  KELLEY 


GOOOARO  SPACE  FLIGHT  CENTER 
GREENBELT,  MARYLAM)  20771 


U.  S.  ARMY  A8ER0EEN  RESEARCH  AM >  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
AfiEROEEN,  MARYLAM) 


OR.  S.  CHANDRA 
OR.  K.  MAE  DO 


OR.  J.  HEIFCRL 
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OR.  M.  BIOAOI 
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CAMERON  STATION 
ALEXANDRIA,  VA.  22314 
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OTHERWISE  2  COPIES)  12CY  ATTN  TC 

UNIVERSITY  OF  CALIFORNIA 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 
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